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Abstract —This paper addresses the problem of identifying a 
linear time-varying (LTV) system characterized by a (possibly 
infinite) discrete set of delays and Doppler shifts. We prove that 
stable identifiability is possible if the upper uniform Beurling 
density of the delay-Doppler support set is strictly smaller than 
1/2 and stable identifiability is impossible for densities strictly 
larger than 1/2. The proof of this density theorem reveals 
an interesting relation between LTV system identification and 
interpolation in the Bargmann-Fock space. Finally, we introduce 
a subspace method for solving the system identification problem 
at hand. 

I. Introduction 

Identification of deterministic linear time-varying (LTV) 
systems has been a topic of long-standing interest, dating back 
to the seminal work by Kailath m and Bello El, and has 
seen significant renewed interest during the past decade na¬ 
il. The problem arises in numerous application areas such 
as radar imaging and wireless communication. The formal 
problem statement is as follows. We want to identify the LTV 
system TL from its response 

Vt G R, (Jlx)(t) := ^ a kx(t — T k )e~ 2m, ' kt (1) 
feen 

to the probing signal x, where (rfc, Vk ) are delay-Doppler shift 
parameters, ak are the corresponding complex weights, and O 
is a possibly infinite index set. 

Contributions: The purpose of this paper is to establish 
fundamental limits on the stable identifiability of "H in 0 
in terms of {ak, rfc, Vk}k^n- Our approach is based on the 
following insight. Defining the complex discrete measure 
Hu := Y^k£n a k$T k ,v k on R 2 , the input-output relation 0 
can be rewritten as 

Vt e R, ( r Hx)(t) = [ x(t — T)e~ 2mvt v). 

Jr 2 

Identifying the system T-L thus amounts to recovering the 
discrete measure fin from ( Hx)(t ), t £ R. This formula¬ 
tion reveals an interesting connection to the super-resolution 
problem as studied by Donoho 0 , where the goal is to 
recover a complex discrete measure on R (a weighted Dirac 
train) from lowpass measurements. The problem at hand, albeit 
formally similar, differs in two important aspects. First, we 
want to recover a measure fiu on R 2 , that is, a measure on 
a two-dimensional set, from observations in one parameter, 
namely (TLx)(t), t £ R. Second, the lowpass observations 
in 0 are replaced by short-time Fourier transform-type ob¬ 
servations, where the probing signal x appears as the window 


function. These differences make for quite different technical 
challenges. Nevertheless, we can follow the spirit of Donoho’s 
work 0, who established necessary and sufficient conditions 
for unique recovery in the super-resolution problem. These 
conditions are expressed in terms of the uniform Beurling 
densities of the measure’s support set and are derived using 
density theorems for interpolation in the Bernstein and Paley- 
Wiener spaces 0 and for balayage of Fourier-Stieltjes trans¬ 
forms S3- Here, we will use a density theorem for interpolation 
in the Bargmann-Fock space m-m. Our main result says 
that stable identifiability is possible if the upper uniform 
Beurling density of the set {(t^, fk)}kGn is strictly smaller 
than 1/2 and stable identifiability is impossible for densities 
strictly larger than 1/2. Moreover, we present a subspace 
method for recovering the parameters {ak, Tk, Vk}k£n from 
Hx when a; is a Gaussian signal. 

Relation to prior work: Kozek and Pfander showed in 0 
that Gabor frame operators on rectangular lattices a -1 Z x 
6 _1 Z are identifiable if and only if ab ^ 1. This prob¬ 
lem is equivalent to the LTV system identification problem 
considered here for the pairs ( Tk. Vk) lying on the lattice 
a _1 Z x 6 _1 Z. The identifiability condition in 0 holds ir¬ 
respectively of whether the set {{Tk,fk)}k £o occupies the 
entire lattice or just parts of it, i.e., the result does not take 
into account the density of {{Tk,Vk)}ken across R 2 . The 
results in 0 pertain to the identification of LTV systems with 
compactly supported spreading function and do not encompass 
operators defined by ©. in mi, a necessary condition for 
identifiability of a set of Hilbert-Schmidt operators defined by 
atomic decompositions is given; the condition is expressed in 
terms of a “2-dimensional” Beurling density, but the operator 
class considered does not contain the operators characterized 
by 0. In 0, it is shown that parametric underspread LTV 
systems, i.e., LTV systems with delay-Doppler spread product 
smaller than 1, can be identified provided that the time- 
bandwidth product of the probing signal is large enough. 
In 0, a convex programming algorithm for stable recovery 
of the triplets {ak, Tk, fk}k^o. from a finite number of noisy 
samples of 'Hx is proposed; the algorithm assumes that the 
points {{rk, izfc)}fcefi obey a minimum separation condition. 

Notation: The complex conjugate of z £ C is denoted by 1. 
For a Hilbert space H, we write (•, ■) H and ||-|| ff for the inner 
product and norm on H, respectively. Linear operators are 
denoted by uppercase calligraphic letters. For v £ R, we define 
the modulation operator {M. v x){t) := a;(f)e -2 '"'' t , and for 


ret, the translation operator {J~ T x)(t) := x(t — r). For a 
measure space (X, E ,p) and a measurable function 93 : X —> 
C, we write f x ip(x)p(x) for the integral of tp with respect to 
p, and we set dx := X(x) if A is the Lebesgue measure. If X is 
a topological space, B(X) designates the Borel cr-algebra over 
X and Xi(X) is the space of all complex Radon measures 
on {X,B{X)). For x £ X, S x £ Xi(X) denotes the Dirac 
measure at x, which for B £ B{X) is given by 5 X (B) = 1, 
if x £ B, and S X (B) = 0, else. The support supp(^t) of a 
complex Radon measure p £ XI (X) is the largest closed set 
CCA' such that for every open set B £ B (X) satisfying 
B fl C f 0, it holds that p(B (T C) 0. We define the total 
variation (TV) norm ||-|| TV on M(X) as ||/r|| TV := |/u|(A), 
where \p\ £ Xi{X) is the total variation of p given by 

VB 6 B( X), \p\(B):= sup £ \p(A)\ 

’ ren ( B ) Aeir 

with 11(73) denoting the set of all partitions of B. For a 
discrete measure p = a *$x k £ M(X), where {x k }k&n 

and {ctfcjfcen are sequences in X and C, respectively, we 
have ||/i|| TV = n |a!fc|, and we define the norm ||/r|| 2 := 

\fj2ken l a fc| 2 - <50 (R) stands for the Feichtinger algebra {15]. 
A set A C R 2 is said to be discrete if for all A £ A 
one can find 6 > 0 such that ||A— A'||^ > 6 for all 

A' £ A\ {A}. A set A C l 2 is said to be uniformly discrete 
if inf{11A - A '|| £2 :A,A'eA,A/ A'} > 0. 

II. Problem formulation 

Throughout the paper, we let S C /1Z 2 be a lattice, where 
A £ l 2x2 is an invertible matrix. We consider LTV systems 
characterized by the bounded linear operator H : So (R) —> 
S’o(R) defined as 

Vi £ R, (Hx)(t) := f x(t — T)e~ 2mvt hu{t, v) , 

J R 2 

where p-y £ XI (E). The vector space of all such operators is 
denoted by Jff(E) and is equipped with the norm \\H\\^, '■= 
||/i^|| 2 . We pursue two principal goals. First, we want to 
establish conditions under which H £ Jf?(E) is identifiable, 
that is, one can find a signal x £ S'o(R), henceforth called 
probing signal, such that Hx uniquely determines 'H. The 
second goal is to find a method for recovering the triplets 
{flfcCfc, Vk}ken associated with H from the reponse of H to 
a probing signal. Throughout the paper, we consider stable 
identifiability (hereafter simply referred to as identifiability) 
which guarantees that TL can not only be recovered from 
Hx, but small errors in Hx also result in small errors in the 
identified operator. 

Definition 1 (Stable identifiability). Let J 1 C J^(S) be a 
set of operators. We say that J 1 is identifiable if there exist 
x £ S 0 (R) and constants C \, C 2 , 0 < C\ ^ C 2 < 00 , such 
that 

Ci ||H - 1CWjg, < ||Hx - lCx \\ L 2 < C 2 ||H - 1CWjg, (2) 

for all H,YC £ J 1 . 


The upper bound in (0 is met trivially thanks to the 
following result. 

Proposition 1. Let x £ <So(R). For all H,1C £ it holds 

that 

\\Hx-1Cx\\ l2 ^\\H-1C\\^\\x\\ l2 . 

Proof: Follows directly by application of [ 16. Thm. 12] 
as x £ So(R) an d supp(py — p/c) C H is uniformly discrete. 


Proving that a set of operators J? C is identifiable 

therefore amounts to finding a probing signal x £ So (R) such 
that the lower bound in 0 holds for all H,1C £ J?. 


III. Main result 

We derive identifiability conditions for operators H £ 
in terms of the uniform Beurling density of 
supp(/u-h) = {(Tfc,u fc )} fce fi. The uniform Beurling density 
measures the average number of points in supp(py) per unit 
cell of R 2 . We will actually have to work with lower and upper 
uniform Beurling densities. 


Definition 2 (Lower and upper uniform Beurling densities, 
0 p. 346] H3 p. 47]). Let X be a uniformly discrete set 
in R 2 . Let Q C R 2 be a compact set of measure A (Q) = 1 
whose boundary has measure 0. For r > 0, let n~{ A, rQ ) and 
n + (A,rQ) be, respectively, the smallest and largest number 
of points of A contained in any translate of rQ. The quantities 


lim inf ■ 

r—foo 


'(A ,rQ) 


and lim sup ■ 


"(A, rQ) 


do not depend on Q and define the lower and upper uniform 
Beurling densities of A, respectively. These quantities are 
denoted by D~ (A) and D + ( A), respectively. If D~{ A) = 
D + { A), then A is said to have uniform Beurling density 
D(A) := D~(A) = D+{A). 


In the remainder of the paper, we will often deal with sets of 
complex numbers {z k }ken, whose lower and upper uniform 
Beurling densities we define to be the lower and upper uniform 
Beurling densities of the sets {(Re { z k }, Im {zfc})}*. e n in R 2 . 
We are now ready to state our main result. 


Theorem 1 (Density criteria for identifiability). Let a be a 
positive number and assume that 5 has uniform Beurling 
density D(E) = l/det(A) ^ 2a. Define the set of opera¬ 
tors j4f a (E) := {H £ Jf’(E): D + (supp(py)) ^ a} and let 
g(t) := v /r Be _7rB t 2 D, ^ g ]R ; where B > 0. 

a) If a < 1/2, then Jlf a (E) is identifiable by g. 

b) If a > 1/2, then J^ a (E) is not identifiable. 


Theorem Q] says that operators H £ J(?(E) are identifiable 
if they are “sparse enough” in the sense of the upper uniform 
Beurling density of supp {pu) satisfying D + (supp(py)) < 
1/2. The proof of Theorem Q] is given in Section ITVl and relies 
on the theory of interpolation in the Bargmann-Fock space 
H(C) of entire functions ip for which 

IMI^:= \<p{z)f e-^dzj V2 


< 00 . 





This powerful theory was developed by Seip and Brekke 
in m-m and led to strong results in Gabor theory es¬ 
tablishing density criteria on sets A := {(Tk,Vk)}k £n for 
{M Vk %- k g}k£Q. to form a Riesz sequence in L 2 (R). 

The result that comes closest to our Theorem Q] is 0 
Thm. 4.1] summarized next. Consider the set of Gabor frame 
operators SP a , b '■= : g £ L 2 (K), £ VF(R)|, where 

a, b > 0, kF(R) is the Wiener space 

:= | / £ L 2 (R): ^ ||/- l[ fe , fc+ i )|| L<50 < oo 1 , 

l fcez J 

and 5“’^ is the linear operator mapping / £ L 2 (R) to 

s g’hf '■= ( ab ) _1 ^2 ( h ^ M m/aTn/ b g) M m/a T n/b f 
ra,n(=Z 

in L 2 (R). y ayb is shown in 0 to be identifiable if and only 
if ab ^ 1. This statement is a universal result in the sense of 
holding irrespectively of whether the set of delay-Doppler shift 
pairs {(t?;, t/fc)}fc S n occupies the entire lattice a -1 Z x & -1 Z 
or just parts of it. In contrast, our result takes into account 
the lattice occupation density. The factor of two difference 
between 0 Thm. 4.1] and Theorem [T| in the critical density 
stems from the fact that for "H,/C £ o5^(El), we have % — K, £ 
M 2 q.(E), whereas the set of operators y a ,b is linear, implying 
that for K,IC £ y a ,b, H — K, £ y a , b as well. 

IV. Proof of Theorem 0 

We first collect some basic facts about interpolation in the 
Bargmann-Fock space m-m. 

A. Preparatory material 

Definition 3 (Set of interpolation, ifTOl l. Let T := {^fejfcgn Q 
C be a discrete set. If for every sequence {wk}ke fi £ m 
there exists a function ip £ -F(C) such that e ^ 2 p(zk) = 
Wk, for all k £ Cl, then T is said to be a set of interpolation 
for B(C). 

Theorem 2 (Density theorem for interpolation in the 
Bargmann-Fock space, lUOl Thm. 1.2], 112] Thm. 1.2]). A 
discrete set T C C is a set of interpolation for T(fZ) if and 
only if it is uniformly discrete and D + (r) < 1 . 

The Bargmann-Fock space is isomorphic to L 2 (R), since 
the Bargmann transform defined for / £ L 2 ( R) as 

/ OO 

f(u)e 2 ' KUZ 

-OO 

is an isometric isomorphism from L 2 (R) to -F(C). The 
Bargmann transform is closely related to the short-time Fourier 
transform (STFT) with respect to a Gaussian window 00 
Def. 3.4.1]. Specifically, the STFT of a; £ L 2 { R) with 
Gaussian window function g[t) := \fBe~ nB 4 / 2 , t £ R, 
defined as 

/ OO 

x(t)g(t — T)e~ 2nlvt At 

-OO 

can be expressed in terms of the Bargmann transform as 

V(r, u) £ R 2 , (V g x)(r, v ) = (Bf)(z)e-^ 2 ' 2 e~™\ (3) 


dw 


where we set f(u) ■= (\/2/BY^xfus/ 1 !/B) and z := 
tB/s/2 — ivs/2/B. 

Definition 4 (Riesz sequence, lfT 8 l Chap. 3]). Let {f n }ne z be 
a sequence of functions in L 2 (R). The sequence {f n }ne z is 
a Riesz sequence in L 2 (M) if and only if there exist constants 
C i, C 2 , 0 < Ci < C 2 < 00 , such that 


/El c -l 2 < 

^ ^ c nfn 


1 n£Z 

nEZ 

L 2 V 


nEZ 

for all {c n } n£ z £ £ 2 (Z). 

An important consequence of the relation between the 
STFT and the Bargmann transform is the following: the set 
{JClv k T Tk g}ken forms a Riesz sequence in L 2 (R) if and only 
if the set T := {z k }k&n, where z k := TkB/y/2- iv k V^/B, 
for all k £ Cl, is a set of interpolation for .F(C). 

B. Proof of Theorem]!} Statement [ajl 

Our proof is inspired by the technique used in 0 Thm. 1.1] 
to prove that discrete complex measures of the form p = 
1 a k<ikT are uniquely characterized by their Fourier trans¬ 
form p (/) = J2kez a k e ~ 2mkTf for I/I < fc if the support of 
p obeys 0 + (supp(^i)) < f c . 

Let a < 1/2 and U, K. £ Jf a (E). If U = 1C, then © holds 
trivially. We consider the case H 7 ^ /C in the following. The 
support of the measure p := pu — pic is contained in A := 
supp(/i-H)u supp(/r;c). Since A C Sis uniformly discrete, we 
have 

D + { A) < D + (supp(p n )) + D + (snpp(pic)) < 2 a < 1 . (4) 

We write A = {(t*,, Ok)}ke n and define the corresponding set 
r := {z k }ken, where z k ■= x k + iyk, x k ■= T k B/\/2, and 
Uk := —OkVZ/B, for all k £ Cl. Let Q := [0, l ] 2 and Qb ■= 
[0, y[2/B\ x [0, B/y/2] and note that A (Q) = X(Qb) = 1 and 
the boundaries of Q and Qb both have measure 0. Moreover, 
for r > 0, (. x k ,yk ) € rQ if and only if (r k ,v k ) £ rQ B - The 
largest number of points of T found in any translate of rQ then 
equals the largest number of points of A found in any translate 
of rQs for r > 0. Therefore, we have D + (T) = D + { A). 
From 01 , it then follows by application of Theorem 0 that 
there exists a function p £ T(C) solving the interpolation 
problem e _ 7 r l 2fe l ^ 2 p(zQ = Wk, for all k £ Cl, where the 
sequence {w k }ken £ £ 2 (Cl) is defined as 

Oik 


\/k £ Cl, Wk = 


_g7T2Tfco I'/cq 


(5) 


\Jj 2 ken \ a k 

with a k := 77 ({(Tfc, )}), for all k £ Cl. Note that ak Y 

0 for at least one k £ Cl, say ko £ Cl, since M / 1C. It 
therefore follows that fZ k . e n \ a k\ 2 > 0 , which ensures that 
the Wk in 0 are well-defined. By ]2ken \ w k\ 2 = 1 < 00 , the 
sequence {wk}ke o is in £ 2 {Ci). Since the Bargmann transform 
B: L 2 (R) —7 J-(C) is an isomorphism and p £ ^(C), there 
exists an / £ L 2 (R) such that <p = Bf. Now, defining x(t) := 
(B/V2y/ 2 f(tB/V2), t £ R, we make use of the relation 0 
between the STFT and the Bargmann transform to write 

Vfc £ 12, p{z k ) = (Vgx)(T k ,u k )e” M 2 /V-^. 










It therefore follows that 


vfc g n, (v g x)(Tk, Vk) = 


By construction of x, it holds that 


Oik 


\Jj2k&n \ a k\ 


(6) 


/ (V s x)(r,^)r?(r,z/) = V a k (V g x)(T k , v k ) = /V |a fc | 2 
•^ R2 fcen y fceo 

= IMI 2 = II Pu - M/c|| 2 = 11^ - ^lljr > 

where the last equality is by the definition of ||-||On the 
other hand, we have 


/R 2 


(V g a:)(T, v)t)(t, v) 


f 


a;(f) fl (f-r)e- 2 ^df U<r,i/> 

-OO / 

[°° x{f) ( f g{t — T)e- Mvt V (T ,«/)) d t (7) 

J— oo \«/R 2 J 

- (/C£f)(f))df 

<||z|| L2 ||77. g -/Cg|| i2) (8) 

where we used Fubini’s theorem to get 0. The conditions for 
Fubini’s theorem are satisfied as 

[ ( / \x(t)g(t-T)e~ 2mi,t \dt\ \v\{t,u) 

«^R 2 V J — OO / 

^ II x IIl 2 llsllz/2 ||g|| T y < oo. 

Fi nally, ([Sj> foll ows from the Cauchy-Schwarz inequality. Since 

\Jj2ken\ a kf < Efcen Mj^JMItv and II^IItv < oo by 

7} £ Af(S), we have ^/Efcen l a fc| 2 < 00 ■ Combining this 
with afc 0 ^ 0, it follows from © that (V 9 x)(rfe 0 , v ko ) ^ 0. 
Since 

V(r, i/) £ R 2 , (V s x)(t, i/) = ^(^)e _7r|2|2/2 e _7riTI/ , 

where z = rB/y/2 — iv y/2/B, and since ip is an entire 
function, V, r c is continuous. This implies that one can find a 
neighborhood V of ( r ko , v ko ) such that (V g a;)(r, v) ^ 0 for all 
{t, v) £ V. Therefore, ||V g x|| 2 2 ^ f v |(V g x)(r, v)\ 2 drdv > 
0. Since ||g|| i2 = 1, we have ||a;|| L2 = l|V 9 a;|| i2 > 0. This 
allows us to write 0 in the form 

Ci\\H-K\\ x ^\\Hg-Kg\\ L *, 
where C\ := 1/ ||a;|| i2 > 0, thereby completing the proof. 


C. Proof of Theorem [7] Statement 15]) 

By contraposition, we show that a ^ 1/2 if Jf? a (E) is iden¬ 
tifiable. To this end, let us assume that one can find a probing 
signal x £ ,3’o(R) and constants C\, C 2 , 0 < C\ ^ C 2 < 00 , 
such that (0 holds for all 7-1, JC € We then construct 

sequences {77 n } n gN an d {^nJngN of operators in J^(E) that 
have Ay n := supp(/r?i n ) and A/c n := supp(gjt„) disjoint 
for all n £ N, and show that if 0 is to hold for 7~L = 77 n 
and K = /C n , for n —> 00 , we necessarily have a ^ 1/2. 
Specifically, with d := D(E.) = l/det(A), we construct 
{ A «.} n£ N as follows. By density of Q in R, we can find, 


for every n £ N, a q n £ Q such that 

ala 

-T - - < Qn < -J. (9) 

d n d 

By assumption, d f 2a, and hence 0 ^ q n ^ 1/2 for all 
n ^ \_a/d\ + 1 =: uq. We can therefore write q n = t n /m n , 
where m n , £ n £ N and 0 < i n ^ m n /2 for all n ^ no- Next, 
define the set A 


n 


:= A Ay , where 

Tt-rt ’ 

t n -1 

: = IJ (rri n Z + r ) 2 

r—0 

is a subset of Z 2 with uniform Beurling density D( Ay n ) = 
q n . The set A^ n is then a subset of E with uniform Beurling 
density D(Ay n ) = q ra /det(A) = q n d ^ a for all n ^ n^. 
Moreover, it follows from 0 that I)(Ay n ) converges to a as 
n —t 00 . We next construct the sequence {A^ n } rae N. Again 
invoking the density of Q in R, we can find q' n £ Q such that 
a — 1/n 




a 


( 10 ) 


(1 - q n )d " ^ ^ (1 -q n )d’ 
for all n £ N. As d ^ 2a and q n ^ 1/2, we have 0 < Qn < 1 
for all n ^ [1/ckJ +1 =: n\. We define A^:„ := AA/c n , where 

m n —1 C _1 

Ak„ : = (J (J ( m nm' n Z + m„s + r) 2 

r=£ n s =0 


is a subset of Z 2 \ with uniform Beurling density 

D(A/c n ) = q' n { 1 — q n )d ^ a. The set A^„ is then a subset of 
S \ Ajc n with uniform Beurling density q' n (l — q n )/ det(A) = 
g^(l — q n )d ^ a for all n ^ n\. It follows from (ITUl) that the 
sequence D(A/c n ) converges to a. Now let n ^ max{no, n 1 }. 
As Ay n and Ay n are disjoint and both have uniform Beurling 
density, the set A n := A-y n U Ay n has uniform Beurling 
density D(A n ) = q n d + q' n (l — q n )d. We write A n = 
{(U',iD^fc,n)}fceSt n i • = ^ Qn' (x kl nik' k ,n) £ Ay n f 

and f l Kn := {k £ f l n : (■ T k , n ,v k ,n ) G Ax: n }. Next, let 
{«fe}feen„ be a sequence in £ 1 (f2 n ) and define the measures 

- = ^ ^ dK-n -= _ ^ ] a k^Tk,n,^k,n- 

The corresponding operators 7f n and /C„ are in 
Evaluating 0 with 77 = 77 n and /C = /C„ then yields 


Cr 





Y! a kA / l l , k 

fcen„ 


<C 2 

L 2 



We can therefore conclude that , n 7/ t: „ x}is a Riesz 
sequence in L 2 (R). By application of Theorem lfl6l Thm. 13- 
d)], it is then necessary that D + (A n ) ^ 1, thus implying 
g„c7 4- 7^(1 — q n )d < 1. Taking the limit n —> 00 yields 
2a ^ 1 and completes the proof. 


V. Identification algorithms 

In wireless or radar applications contributions to the re¬ 
ceived signal corresponding to distant reflectors and/or scatter¬ 
ed can often be neglected thanks to path loss. It is therefore 
sensible to take the delays and Doppler shifts to lie on a 
compact set [T m i n ,T max ] x [i/ m in, ^max]■ Combined with the 















standing assumption of the (rk,Vk) being supported on a 
lattice, the input-output relation (QJ reduces to 

K 

Vf e R, (' Hx)(t ) = ^ akx(t — Tfc)e _27I ’ Il/fct 

fe=i 

with K finite. This implies that the set supp (fin) is finite and 
hence Z? + (supp(/u^)) = 0. Theorem Q] then tells us that H 
can be identified by the probing signal 

VteR, x{t) :=VBe X p(- nB2{t ~ T)2 y 

Furthermore, in practice, we only have access to samples of 
Hx restricted to a finite time interval, say [r m ; n , T m ; n +T]. The 
identification problem therefore reduces to the identification 
of TL from the samples r m := (' Hx)(t m ) with t m '■= 
T m i n + mT/M for m £ {0,1 ,,M — 1}. Straightforward 
manipulations reveal that 

K 

Tm — ^M—m ^ ^ ^k^k ) 

k=1 


where we set 


, rs ( 7T B 2 T 2 (M — m) 2 \ 

\ M -rn ■■= Vfiexp^- — - j 


otk '■= a,k exp I - 


(rfc r min) ^ — InivkT2niVkT m 


z k ■= exp - 


■nB"T{rk T m j n ) \ 27 T iv k T/M 


M 


for m £ {0,1,..., M — 1}. Therefore, letting y m := 
r M -m/Kx, we obtain 

K 

Vm=Y, ClkZ k l 

fc= 1 


for all in £ {1,2,..., M }. We have therefore reduced the LTV 
system identification problem to the problem of identifying 
cisoids. A variety of subspace algorithms like ESPRIT fl9l . 
MUSIC J20), or the matrix pencil method ITTl can be used 
to solve this problem, provided that M ^ K + 1. For an in- 
depth discussion of subspace cisoid estimation algorithms we 
refer the reader to li22l . In the presence of noise, the success 
of subspace methods is highly dependent on the condition 
number of the Vandermonde matrix with poles {zk}l 1231 . 
In HQ, it was shown that these Vandermonde matrices are 
relatively well conditioned, provided that the nodes lie inside 
the unit disk, remain close to the unit circle, but are not 
extremely close to each other. 
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